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Abstract. Using Hartree-Fock orbitals with residual Coulomb repulsion, we study spinless fermions in a 
two dimensional random potential. When we increase the system size L at fixed particle density, the size 
dependence of the average inverse compressibility exhibits a smooth crossover from a 1/L 2 towards a 1/L 
decay when the Coulomb energy to Fermi energy ratio r s increases from to 3. In contrast, the distribution 
of the first energy excitation displays a sharp Poisson-Wigner-like transition at r s m 1. 



PACS. 71.30.+h Metal-insulator transitions and other electronic transitions 
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72.15.Rn Quantum local- 



1 Introduction 

The interplay of disorder and interactions in two dimen- 
sional (2c?) electronic systems is currently a central prob- 
lem in condensed matter physics M . The importance of in- 
teractions is illustrated by several phenomena, notably the 
discovery of a metallic phase in various two-dimensional 
devices ||. Since the metal-insulator transition occurs at 
a Coulomb energy to Fermi energy ratio r s w 10, Coulomb 
repulsion should be crucial in understanding this phe- 
nomenon. Compressibility measurements indicate that the 
high r s insulating phase is incompressible |3) and spatially 
inhomogeneous ||. 

In weakly disordered quantum dots, the inverse com- 
pressibility gives the spacing between adjacent conduc- 
tance peaks. The peak spacing statistics, obtained in ex- 
periments in which r s ss 1, displays fluctuations larger 
than those predicted from Random Matrix Theory, char- 
acterized by a Gaussian-like distribution instead of a Wigner- 
Dyson distribution |^|,^U^]. This suggests a breakdown of 
the naive single particle picture, assumed in the constant 
interaction model. 

In strongly disordered insulators, the long range na- 
ture of interactions leads to a soft Coulomb gap in the 
single particle density of states at the Fermi level M , ob- 
served in tunneling spectroscopic experiments in 3c? H and 
2c? |h]] nonmetallic semiconductors. Assuming that single 
electron hoppings dominate the transport, the Coulomb 
gap leads to a crossover in the temperature dependence of 
the resistivity p(T) from the Mott variable range hopping 
law (p(T) = pMexpiTo/T) 1 / 3 in dimension d = 2) to the 
Efros-Shklovskii behavior (p(T) = p E sexp(T ES /T) 1 ' 2 ) 
||. This crossover has been reported jll| not only de- 



creasing temperature but also decreasing carrier density 
(around r s w 1.7). 

In this paper, we numerically investigate the physics 
of 2c? Coulomb interacting spinless fermions in a random 
potential for relatively low values of the factor r s (< 3), 
when the single particle spectra display either Wigner- 
Dyson statistics for weak disorder (diffusion) or Poisson 
statistics for strong disorder (Anderson localization). As 
exact diagonalization studies are restricted to small sys- 
tem sizes [^2), we use a numerical method [ |l3lJl4|| famil- 
iar in quantum chemistry as the configuration interaction 
method (CIM) The method consists in diagonalizing 
the Hamiltonian in an energetically truncated basis built 
of the low-energy states of the corresponding Hartree-Fock 
(HF) Hamiltonian. This method allows us to study the 
ground state and the lowest energy excitations for differ- 
ent system sizes L at constant electronic density. 

We show that the average inverse compressibility A2 
exhibits a smooth crossover from a 1/L 2 towards a 1/L 
decay when r s increases from to 3. At the same time, 
the distribution of A2 evolves from the Wigner-Dyson dis- 
tribution (or the Poisson distribution, if disorder is strong 
enough or L large enough) towards a Gaussian-like shape. 
Therefore fluctuations in A2 are determined by Coulomb 
repulsion rather than from single particle level fluctua- 
tions. 

We have also studied the energy level spacing between 
the ground state and the first excitation for strongly dis- 
ordered clusters: the average indicates a smooth opening 
of a gap when r s increases, while the distribution exhibits 
a sharp Poisson-Wigner-like transition at r s — rf 1=3 1.2. 
The critical threshold is characterized by a scale in- 
variant gap distribution, reminiscent of the one particle 
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problem |16 at a mobility edge. However, it is only the 
distribution of the first spacing which exhibits such a tran- 
sition, the distributions of the next spacings remain Pois- 
sonian and are essentially unchanged when r s varies. 

The paper is organized as follows: in Section we in- 
troduce the model, in Section || we discuss the numerical 
technique used and its limits of validity; in Section |i| we 
present our compressibility data; in Section |H| we discuss 
the statistical properties of the lowest energy excitations; 
in Section ^| we discuss the implications of our results for 
hopping conductivity experiments; in Section ^ we present 
our conclusions. 



2 The model 

We consider a disordered square lattice with M = L 2 sites 
occupied by N spinless fermions. The Hamiltonian reads 



H 



"* £ C * C J + £ ViUi + U Y1 



<hj> 



(1) 



where c\ (q) creates (destroys) an electron in the site i, the 
hopping term t between nearest neighbours characterizes 
the kinetic energy, Vi the site potentials taken at random 
inside the interval [— W/2, +W/2], rij = c\ci is the occu- 
pation number at site i and U measures the strength of 
the Coulomb repulsion. The boundary conditions are pe- 
riodic and Tij is the inter-particle distance for a 2c? torus. 
If a* B = h 2 e/(m*e 2 ), to*, e, a and n s — N/(aL 2 ) denote 
respectively the effective Bohr radius, the effective mass, 
the dielectric constant, the lattice spacing and the carrier 
density, the factor r s is given by: 



since in our units 
n e = N/L 2 . 



U 



h 2 /{2m*a 2 ) -> t, e 2 /(ea) 



(2) 



U and 



3 Configuration interaction method 

A numerical study of the model ([!]) via exact diagonal- 
ization techniques for sparse matrices is possible only for 
small systems and does not allow us to vary L for a 
constant density. We are obliged to look for an approxi- 
mate solution of the problem, using the Hartree-Fock or- 
bitals, and to control the validity of the approximations. 
One starts from the HF Hamiltonian where the two-body 
part is reduced to an effective single particle Hamiltonian 
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Fig. 1. Relative error in the HF ground state energy as a func- 
tion of r s , for N = 4, L = 6, W = 5 (circles), W — 8 (squares), 
and W — 15 (diamonds). The contribution of a positive uni- 
form background has been included in the energies. Brackets 
indicate disorder average (here over 10 2 configurations). 



system sizes the single-particle problem (||) is still non- 
trivial, since the self-consistent iteration can be trapped 
in metastable states. This limits our study to small r s and 
forbids us to study by this method charge crystallization 
discussed in j[2) at a larger ~ 10 . 

Fig. |l| shows that the HF approximation gives a good 
estimate of the ground state energy. This approximation 
becomes worse for (1) larger interactions or (2) smaller 
disorder. The first characteristic is due to the fact that 
the ground state HF energy is exact in the limits U — * 
and U — > oo, in which the ground state is given by a sin- 
gle Slater determinant, but deviations from this simple 
picture are expected at intermediate U values. The lat- 
ter feature can be understood as complicated many-body 
screening effects, which are effective up to a distance of 
the order of the localization length, cannot be reproduced 
within simple HF. 

The main advantage of the HF approximation is that 
it reproduces well the single particle density of states |L9|, 
particularly the soft Coulomb gap at the Fermi energy j. 
However, the approximations involved in the HF method 
are uncontrolled. The mean field HF results can be im- 
proved using the configuration interaction method fl^ , |l4] , |i"5"f 
Once a complete orthonormal basis of HF orbitals has 
been calculated: 



5 r * 



-cla 



(3) 



where (...) stands for the expectation value with respect 
to the HF ground state, which has to be determined self- 
consistently. For large values of the interaction and large 



H HF (\ipi 



(4) 



with a = 1, 2, . . . , L 2 , it is possible to build up a Slater 
determinants' basis for the many-body problem which can 
be truncated to the Njj first Slater determinants, ordered 
by increasing energies. The two-body Hamiltonian can be 
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Fig. 2. CIM result (filled circles) compared to HF approxi- 
mation (empty circles) for N = 9, L = 12, W = 15. Mean gap 
< AEq > and gap distribution at r s = 4.5 (inset). Disorder 
average is over 10 4 configurations. 

written as 

H int = \ ^2 Qlp^a d l d sdy, (5) 

with 

tj \ ^ 

0£^3. ^^^^ _ m 

and S a = J^j ^aif) ]^)- O ne gets the residual interaction 
subtracting Eq. || from Eq. |[ This keeps the two-body 
nature of the Coulomb interaction, and if TV ^ 2 it is still 
possible to take advantage of the sparsity of the matrix 
and to diagonalize it via the Lanczos algorithm. 

Fig. U compares HF and CIM results. Labelling the 
iV-body levels (i — 0, 1, 2, ...) by increasing energy we 
have studied the first spacing AEq = E\ — Eq. When 
r s > 1, the residual interaction significantly reduces the 
mean gap, and slightly changes its distribution. Within 
HF approximation, the first excitation is a particle-hole 
excitation starting from the ground state: The mean gap 
reduction is due to correlation effects beyond the particle- 
hole interaction. For example, the electron (hole) polar- 
izes its neighborhood as it creates fluctuations in the local 
charge density. 

We have checked that CIM results agree with the re- 
sults given from exact diagonalization with an accuracy of 
the order 2% when one takes into account the Nh = 10 3 
lowest energy Slater determinants when N = 4, L = 8, 
W = 15, and r s = 5. This means that a basis spanning 
only 0.2% of the total Hilbert space is sufficient for study- 
ing the first excitations. For larger L, exact diagonaliza- 
tion is no longer possible, but one can look the variation 



Fig. 3. First (circles), second (squares), third (diamonds), and 
fourth (triangles) average spacing as a function of the size 
Nh of the truncated Hilbert space, for N — 16, L = 16, 
W = 15, r s — 2.8. The first gap is fitted by the curve 
< AE >= 0.195 + 0.13iV^ - 26 ; therefore one can estimate 
the relative error at N H = 2 x 10 3 as < 10%. Disorder average 
is over 10 2 configurations. 

of the results when Nh increases. In the worst case con- 
sidered (N — 16, L — 16, W = 15, r s — 2.8) the accuracy 
in the first four spacings can be estimated as better than 
10% when N H = 2 x 10 3 (see Fig. §. 

Therefore the CIM method allows to study low energy 
level statistics for r s < 3. However, its accuracy is not suf- 
ficient to determine more sensitive quantities, like a small 
change of the ground state energy when the boundary 
conditions are twisted (i.e. the persistent currents). 



4 Electronic compressibilty 

In quantum dots experiments in the Coulomb blockade 
regime the spacing between consecutive conductance peaks 
is given by 

A 2 (N) = Eq(N + 1) - 2E (N) + Eq(N - 1), (7) 

with E (N) TV-body ground state. This quantity is the 
discretized second derivative of the ground state energy 
with respect to the number of particles, i.e. the inverse 
compressibility. In the constant interaction model, which 
ignores fluctuations in the Coulomb energy 

s e 2 

Qlf3 = -Q S ocyS 0S , (8) 

where C is the capacitance of the dot. This gives Eq(N) = 

e- 2 N(N—l)/2C+J2k=i Vk, where % is the k-th single par- 
ticle energy at U = 0. Substituting this expression into Eq. 
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(|) one finds A 2 = e 2 /C + Vn+i — Vn- Random Matrix 
Theory describes the single particle level spacings when 
the disorder is weak. The peak spacing distribution is ex- 
pected to follow the Wigner-Dyson distribution with fluc- 
tuations in A% whose standard deviation is given by bA 2 = 
y/< A{ > - < A 2 > 2 = - 1 < A >« 0.52 < A >, 

where < A > is the single particle mean level spacing. 
However, experiments [SpLItI performed at r 8 « 1 found a 
distribution which is Gaussian-like and has a larger width, 
up to 7.5Z\ 0. This suggests that fluctuations in A 2 are 
dominated by electron-electron interactions. Several ob- 
served features of the peak spacing distributions have been 
reproduced using exact numerical diagonalization of the 
Hamiltonian model (Q) [§, HF calculations and 
a random matrix model for interacting fermionic systems 
II- 

In this paper we study the inverse compressibility at 
different system sizes, for a constant filling factor n e = 
1/9. We consider N = 4, 9, 16 particles on square lattices 
of size L = 6, 9, 12 respectively. A Fermi golden rule ap- 
proximation for the elastic scattering time gives pi]], for 
n e <C 1, kpl w 192irn e (t/W) 2 , kp and I denoting the 
Fermi wave vector and the elastic mean free path, re- 
spectively. We consider W — 5, 8, 15, corresponding to 
kpl ~ 2.7, 1,0.3; the first case corresponds to a diffusive 
system (L > //aw 2.3, with a lattice spacing), the last 
one to a strongly localized system. 

The following inverse compressibility data are obtained 
with the CIM. We have checked that the residual inter- 
action does not change qualitatively HF results. This is 
consistent with Fig. fit HF approximation gives a good 
estimate of the ground state energy and the inverse com- 
pressibility is a physical observable which depends only on 
the ground state energies at different number of particles. 
Neither a precise knowledge of the ground state wavefuc- 
tion (as in the calculation of persistent currents) nor ex- 
cited states energies (as in studies of spectral statistics) 
are required. 

Fig. [l| shows the L-dependence of the average inverse 
compressibility, which is well fitted by the power law < 
A 2 (r s ) >oc L~ a ( Ts \ with a(r s ) going from 2 to 1 when 
r s goes from to 3 approximately. The value a = 2 is 
expected without interaction (< A% >—< A >oc 1/L 2 ). 
The exponent a — 1 can be understood as follows: for r s > 
1, A 2 is dominated by the charging energy and therefore 
< A% >oc 1/C cx 1/L. In the strongly localized regime 
(W = 15) the 1/L law is justified due to the Coulomb gap 
in the single particle density of states ^]: According to 
Koopmans' theorem ^3] , assuming that the other charges 
are not reorganized by the addition of an extra charge, 

A 2 w e N+1 - e N cx — (9) 

due to the Coulomb gap, where eu is the energy of the 
fc-th HF orbital at a fixed number N of particles. 

Fig. U shows the distribution of inverse compressibili- 
ties, for L = 12, W = 5 and W = 15 (inset). At r s = 0, 
A 2 distributions are close to the Wigner surmise Pw(s = 
Ail < A 2 >) = (7rs/2)exp(-7rs 2 /4) for W = 5 and to 




Fig. 4. Inset: size dependence of the average inverse com- 
pressibility < Ai >, for W — 15, filling factor n e = 1/9, 
r 8 — (circels), 0.4 (squares), 0.8 (diamonds), 1.7 (triangles 
up), and 3.4 (triangles down). Straight lines are power law 
fits < A 2 (r s ) >cx L- a( - r ° ) . Main fi gure: Exponent a(r s ) for 
n e = 1/9, W = 5 (circles), W — 8 (squares), and W — 15 (di- 
amonds). Disorder average for the compressibility data is over 
10 3 configurations. 

the Poisson distribution Pp(s) = exp(— s) for W = 15. In 
the latter case, deviations from the Poisson distribution 
at small s values arc due to the finite system size. In both 
cases at r s = 2.5 distributions show a Gaussian shape. 
This can be understood within the Koopmans' theorem, 
as the HF energies are given by 

N 

e k = {MHMk) + £ (Qak - Qll) , (io) 

a=l 

with Hi one-body part of the Hamiltonian (|l|): Due to the 
small correlations of eigenf unctions in a random potential, 
one can reasonably invoke the central limit theorem (see 
Ref. |22| for a more detailed discussion). 

Fig. [| shows the r s -dependence of the inverse com- 
pressibility standard deviation SA 2 , normalized to the sin- 
gle particle mean level spacing < A >. One can see that 
5A 2 / < A > increases slowly for r s < 1 and then more 
rapidly. This results are in contrast with random phase 
approximation estimates |25[|2r| , giving fluctuations of the 
order A. On the contrary, we have checked that 8A 2 / < 
A 2 > has a weak r s dependence for r s > 1, in agreement 
with findings of Ref. H . This confirms that peak spacing 
fluctuations are dominated by capacitance fluctuations in- 
stead of single particle fluctuations. 

We point out that our data are limited to spinless 
fcrmions. Interesting interaction-induced spin effects have 
been recently reported in experiments |g7| and theoreti- 
cally investigated |||,^,^o| . 
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Fig. 5. Distribution of the normalized inverse compressibili- 
ties for N — 16, L — 12, W — 5, r s = (empty circles) and 
r s = 2.5 (filled circles), fitted by a Gaussian of standard devia- 
tion a = 0.18 (dotted line). Solid and dashed lines give Wigner 
and Poisson distribution respectively. Inset: same at W = 15; 
the Gaussian fit gives a — 0.30. 



Fig. 7. Size dependence of the average gap (first spacing < 
AEo ><x L-" (rs) ), for W = 15, filling factor n e = 1/16. From 
bottom to top: r s = 0, 0.6, 0.9, 1.1, 1.4, 1.7, 2.3, 2.8. Inset: p(r„) 
(circles, characterizing < AEo > and squares, characterizing 
< AEi >, with an average over i = 1 — 3). Disorder average is 
over 10 4 configurations. 




Fig. 6. Inverse compressibility standard deviation SA2, nor- 
malized to the single particle mean level spacing < A > as a 
function of r s , for W = 5, n e = 1/9, L = 6 (circles), L = 9 
(squares), and L — 12 (diamonds). Inset: same at W — 15. 



5 First energy excitation 

We have calculated the first TV-body energy levels Ei, 
(i = 0,1,2,...) for different sizes L, with a large disor- 
der to hopping ratio W/t =15 imposed to have Anderson 
localization and Poissonian spectral statistics for the one 
particle levels at r s = when L > 8. We studied N = 4, 9, 
and 16 particles inside clusters of size L = 8, 12, and 16 
respectively. This corresponds to a constant low filling fac- 
tor n e — 1/16. We studied an ensemble of 10 4 disorder 
configurations. 

The first average spacing < AEq > is given in Fig. [?|. 
It exhibits a power law decay as L increases, with an ex- 
ponent [3 given in the inset. One finds for the first spacing 
that (3 linearly decreases from d = 2 to 1 when r s in- 
creases from to 3. The next mean spacings < AEi >=< 
Ei + \ — Ei > depend more weakly on r s , as shown in Fig. 

For r s — 0, the distribution of the first spacing s = 
AEq I < AEq > becomes closer and closer to the Pois- 
son distribution -Pp(s) when L increases, as it should be 
for an Anderson insulator. For a larger r s , the distribu- 
tion seems to become close to the Wigner surmise Pw{s) 
characteristic of level repulsion in Random Matrix The- 
ory, as shown in Fig. [s] for r s — 2.8 and L — 16. To study 
how this P(s) goes from Poisson to a Wigner-like distribu- 
tion when r s increases, we have calculated a parameter 77 
which decreases from 1 to when P(s) goes from Poisson 
to Wigner: 



var(P(s)) — va,r(P w (s)) 
var(Pp(s)) — v&r(P w (s)) ' 



(11) 
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Fig. 8. Gap distribution P(s) for r 3 — (empty circles) and 
r s = 2.8 (filled circles) when N = 16, L = 16, W = 15, 
compared to Pp{s) and Pw(s). Inset: size invariant P(s) at 
ta 1.2 ; L — 8 (circles), 12 (squares), and 16 (diamonds). 



where var (P(s)) denotes the variance of P(s), var(Pp(s)) = 
1 and vax(Pw(s)) = 0.273. In Fig. |9[ one can see that 
the three curves rj(r s ) characterizing the first spacing for 
L = 8,12,16 intersect at a critical value « 1.2. Our 
data suggest that for r s < the distribution tends to 
Poisson in the thermodynamic limit, while for r s > 
it tends to a Wigner-like behavior pi| . At the threshold 
rf , there is a size-independent intermediate distribution 
shown in the inset of Fig. |[ exhibiting level repulsion 
at small s followed by a exp(— as) decay at large s with 
a w 1.52. This Poisson- Wigner transition characterizes 
only the first spacing, the distributions of the next spac- 
ings being quite different. The inset of Fig. || does not 
show an intersection for the parameter rj calculated with 
the second spacing. The second excitation is less localized 
than the first one when r s = 0, since the one particle lo- 
calization length weakly increases with energy. It is only 
for L = 16 that the distribution of the second spacing 
becomes close to Poisson without interaction, and a weak 
level repulsion occurs as r s increases. 

The observed transition, and the difference between 
the first spacing and the following ones is mainly an effect 
of the HF mean field. For the first spacing, the curves 77 
calculated with the HF data are qualitatively the same. 
At the mean field level the low energy excitations are 
particle-hole excitations starting from the ground state. 
The energy spacing between the first and the second ex- 
cited states is given by the difference of two particle-hole 
excitations and a Poisson distribution follows if the low 
energy particle-hole excitations are uncorrelated. 

Notice that the energy of an electron-hole pair is given 
by ej — ej — U /r»j and the classical argument for the exis- 
tence of a gap in the single particle density of states does 



Fig. 9. Parameter T](r 3 ) corresponding to the first spacing 
AEo at n e — 1/16, W = 15, L — 8 (circles), 12 (squares), and 
16 (diamonds). Inset: rj(r s ) for the second spacing AE\. 



not apply p|. For example, in the constant interaction 
model one pays a gap e 2 /C in the single particle density 
of states but not in the many-body excitation spectrum 
(AE = AE (U = 0) cx 1/L 2 ). Therefore the observed 
opening of a gap for the first energy excitation is a remark- 
able phenomenon beyond the predictions of the classical 
Coulomb gap model. 



6 Hopping conductivity 

In disordered insulators, a crossover || in the temper- 
ature dependence of the resistivity p(T) is induced by 
Coulomb repulsion from the Mott variable range hop- 
ping law (p(T) = pM exp(Xb/T) 1//3 in 2d) to the Efros- 
Shklovskii behavior \p(T) = p E s exp(T E s/P) 1/2 )- 

The usual argument is to consider the length L(T) 
which maximizes the probability P^ of a single particle 
hop from a state with energy to a state with energy ej : 



Pij oc exp 



2L 
T 



kT 



(12) 



with £ localization length. 

To measure possible delocalization effects, we have cal- 
culated for W = 15 the number of occupied sites per par- 
ticle £ s = N/^Ph where pi = (^oKl^o) is the charge 
density of the ground state at the site i. Around r s w 1.2 
and after ensemble average, the maximum increase of £ s 
compared to r s = is negligibly small (2%). Therefore 
noticeable effects come mainly from the mean and the 
distribution of the single particle density of states around 
the Fermi level. 

At low temperatures, hopping conductivity involves 
states near the Fermi level. If one takes for the Coulomb 
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gap its average value £n+i — ejv ~ A% ~ (A + Br s )/ L at - r "^ 
(see Fig. |j), one obtains for the hopping resistivity a smooth distribution at 
and continuous crossover from Mott to Efros-Shklovskii 
hopping, given by: 



interaction-induced Poisson-Wigner-like transition in its 



1.2, with a scale invariant distribu- 



where 



p(T) oc exp 



T(r s ) 



T(r s ) 
T 



l/(a+l) 



A + Br s 



(13) 



(14) 



This prediction neglects the crossover from the Poisson 
distribution to a Gaussian distribution in the inverse com- 
pressibility (see Fig. |E|), which could be included by con- 
sidering for A 2 a more typical value than its average, for 

instance the value A2Q)) for which J^ 2 ^ P(A 2 )dA 2 = 
b, with 6=1/2 for example. This would introduce a 
crossover in T(r s ) around r,«l. 

The existence of a critical r s value for the opening 
of the Coulomb gap can be understood similarly to J|2| . 
The single particle density of states around the Fermi en- 
ergy E F is given by g p{E) rs \E - E F \/U 2 and the 
gap size A g = \E g — Ep\ can be estimated from the 
condition p{E g ) rs ~p, with ~p rs 1/W mean density of 
states for W >■ t, obtaining A g rs U 2 /W. According 
to Fermi golden rule, the inverse lifetime of a Slater de- 
terminant built from electrons localized at given sites is 
r t rs t 2 {l/W){N/L 2 ), with N/{WL 2 ) density of states di- 
rectly coupled by the hopping term of the Hamiltonian 
(|l|). This argument suggests that quantum fluctuations 
can melt the Coulomb gap for r t rs A g , giving r s rs 1. 
Therefore a crossover from Efros-Shklovskii to Mott hop- 
ping conductivity is expected not only increasing temper- 
ature but also increasing carrier density, as observed in 
at r s rs 1.7. 

Notice that in the above arguments we have consid- 
ered only the excitations leading to the Efros-Shklovskii 
behavior for the hopping conductivity, that correspond to 
the withdrawal of one electron or its addition to the sys- 
tem. However, a single electron hop may reorganize the 
location of the other particles, inducing complex many 
particle excitations (see Fig. ||, showing that simple HF 
is unable to reproduce the mean first many-body energy 
spacing when r s > 1). Therefore one cannot exclude that 
correlated hops involving more than one particle affect the 
hopping conductivity m. 



7 Conclusions 

In summary, we have analyzed the inverse compressibil- 
ity and the first energy excitation statistics for spinless 
fermions in weakly and strongly disordered squared lat- 
tices when r s increases from to 3. On one hand, we have 
found a crossover in the electronic inverse compressibil- 
ity from a 1/L 2 towards a l/L decay; at the same time 
its distribution evolves towards a Gaussian shape. On the 
contrary, our data for the first many-body energy spacing 
are consistent with a smooth opening of a gap and a sharp 



tion at the critical point showing intermediate statistics. 
For those small values of r s , the ground state energy is es- 
sentially given by an effective Hartree-Fock mean field. We 
underline that the intermediate quantum regime discussed 
in Ref. jl2| and charge crystallization occur for larger r s 
factors, where the mean-field approximation breaks down 
and the configuration interaction method ceases to ef- 
ficiently converge. This does not allow us to study the 
regime characterizing the experiments reported in Refs. 
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